The Curie temperature and exchange energy between two sublattices in half-metalHc 

greigite Fe3S4 
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High-temperature magnetic measurements have been carried out in hydrothermally synthesized 
greigite (Fe3S4). We show that the Curie temperature of greigite is signihcantly lower than that for 
its iron oxide counterpart Fe304. The lower Tc value (about 677 K) of greigite is in quantitative 
agreement with that calculated using the exchange energy (3.25 meV) and the spin values of the 
two sublattices, which are inferred from the neutron and magnetization data of high-quality pure 
greigite samples. We further show that, with an effective on-site Hubbard energy Ueff = 1.16 eV, 
the lattice constant and two sublattice spins predicted from ab initio density-function theory are in 
nearly perfect agreement with the measured values. The parameter Ueff = 1.16 eV ensures Fe3S4 
to be an excellent half-metallic material for spintronic applications. 



Greigite (Fe3S4) that was first discovered in lake sedi- 
ments from California [l|] is an iron thiospinel. Greigite 
was found to have the inverse spinel structure [ij like its 
iron oxide counterpart, magnetite (Fe304). Greigite has 
been used as a recorder of the ancient geomagnetic field 
and environmental processes and is important for pale- 
omagnetic and environmental magnetic studies ^]^. It 
is also widespread in magnetotactic bacteria that pro- 
duce greigite magnetosomes 0-0] • Therefore, greigite is 
of general interest in geophysics and biology. 

On the other hand, greigite has been less known to 
physicists and material scientists because pure greigite 
samples are difficult to synthesize and some of its funda- 
mental magnetic properties are still unknown. The pre- 
cise value of the Curie temperature Tc remains unknown 
despite variable estimations since 1974. The difficulty 
in determination of the Curie temperature is due to the 
fact that greigite is chemically unstable at high temper- 
atures even in argon environment [8|. High-temperature 
magnetic measurements often revealed chemical decom- 
position that precluded from determination of the Curie 
temperature [1, Q. Spender et al. estimated Tc 

to be 606 K by extrapolating thermomagnetic curves 
to high temperatures. 



11 



made 



Vandenberghe et al. 
Mossbauer spectroscopic measurements up to 480 K and 
extrapolated the effective field of iron in tetrahedral sites 
to obtain a Tc value of at least 800 K. More recent mag- 
netic measurements on high-quality pure greigite samples 
(l^ suggested a Curie temperature of higher than 630 
K. Therefore, the reported Curie temperature of greigite 
ranges from 600 K to 800 K. 



Another important magnetic property is the saturation 
magnetization which is associated with the spin values of 
(Fe'^+)A and (Fe'^+Fe^"'")^ on the respective tetrahedral 
(A) and octahedral (B) sublattices. Coey et al. pjj] sug- 
gested that greigite should have a net magnetic moment 
m > 4: fiB per formula unit considering the spin-only val- 
ues for the ionic moments. They could not account for 
the measured magnetic moment of 2.2 fis- Recent neu- 
tron diffraction experiments on high-quality pure greigite 



samples IJ] indicate a much larger net magnetic mo- 
ment of 3.4 hb- Devey et al. [15[ studied the electronic 
and magnetic behaviors of this material using ab ini- 
tio density-functional theory in the generalized gradient 
approximation (GGA) with the on-site Hubbard energy 
Ueff parameter (GGA -I- U). They calculated the lattice 
constant a and two sublattice moments (m^ and ms) as 
a function of Uef / . Comparing the measured lattice con- 
stant with their calculated results [l^, they found that 
Ueff — 1 eV. This Ueff value leads to rriA niB — il^B, 
in quantitative agreement the measured values from neu- 
tron diffraction [14i|. 

Here we report high-temperature magnetic measure- 
ments on hydrothermally synthesized greigite. We show 
that the Curie temperature of greigite is significantly 
lower than that for its iron oxide counterpart Fe304. The 
lower Tc value (about 677 K) of greigite is in quantita- 
tive agreement with that calculated using the exchange 
energy (3.25 meV) and the spin values of the two sub- 
lattices, which are inferred from the neutron [3] and 
magnetization fli*! data of high-quality pure greigite sam- 
ples. We further show that, with an effective on-site Hub- 
bard energy Ueff = 1.16 eV, the lattice constant and two 
sublattice spins predicted from ah initio density-function 
theory are in quantitative agreement with the mea- 
sured values. The parameter Ueff = 1-16 eV ensures 
Fe3S4 to be an excellent half-metallic material for spin- 
tronic applications. 

Greigite samples were synthesized under a magnetic 
field with a hydrothermal route [l6j . Briefly, ammo- 
nium iron sulfate hexahydrate (0.2 mmol) and cysteine 
(1 mmol) were dissolved in distilled water (25 mL) and 
stirred vigorously for 50 min. The solution was then put 
into a 30 mL Teflon vessel with two circular magnets 
(2.5 kOe) at the top and bottom in a stainless steel auto- 
clave. After 12 h reaction at 180 °C, black products were 
obtained. These products were washed several times in 
distilled water and pure ethanol to remove impurities. 
These samples were dried in vacuum at 60 °C for at least 
6 h. 
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FIG. 1: X-ray difFraction (XRD) spectrum for a greigite sam- 
ple synthesized under a magnetic field of 2.5 kOe. All the 
peaks can be indexed by the spinel Fe3S4 phase with a lattce 
constant a — 9.877 A. The mean diameter of the crystallites 
(subparticles) is determined to be about 36 nm from the width 
of the (311) peak. 



Samples were checked by x-ray difFraction (XRD) right 
after synthesis. Fig. la shows XRD spectrum for a greig- 
ite sample synthesized under a magnetic field of 2.5 kOe. 
All the peaks can be indexed by the spinel Fe3S4 phase. 
The lattce constant a is evaluated to be 9.877 A, in excel- 
lent agreement with the reported value (9.876 A) [H [l^. 
The mean diameter of the crystallites (subparticles) is 
determined to be about 36 nm from the width of the 
(311) peak. 

Magnetization was measured using a Quantum Design 
vibrating sample magnetometer (VSM). The moment 
measurement was carried out after the sample chamber 
reached a high vacuum of better than 9x 10~^ torr. The 
absolute measurement uncertainties in temperature and 
moment are less than 10 K and 1 x 10^^ emu, respectively. 
We used the same heating and cooling rate of 30 K/ min 
for thermomagnetic measurements, which can ensure a 
thermal lag of less than 10 K. 

In Fig. 2, we plot magnetic hysteresis loop at 300 K 
for the greigite sample. Since the magnetization is al- 
most saturated at a magnetic field of 20 kOe, we can 
take the magnetization at 20 kOe to be the saturation 
magnetization Ms- The Ms value is found to be about 
40 emu /g, which is close to that of the similarly prepared 
sample |16| but significantly below the value (59 emu/g) 
reported for high-quality greigite samples with d = 14 ^m 
(Ref. 1121 )• The reduction of the Ms value in our sample 
may be due to both finite-size effects 



17| and possible 



presence of nonmagnetic amorphous phases. 

Figure 3a shows the temperature dependence of the 
magnetization for the greigite sample, which was mea- 
sured in a magnetic field of 10 kOe. The warm-up data 
up to 780 K indicate a magnetic transition around 700 
K. A substantial magnetization between 720 and 780 K 
implies that a small fraction of the Fe3S4 phase was con- 
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FIG. 2: Magnetization versus magnetic field at 300 K. The 
saturation magnetization Ms is found to be about 40 emu/g. 



verted to the magnetic Fe304 phase in this temperature 
regime. Since this sample was measured after the VSM 
high-vacuum cryopump was exposed to air during a ser- 
vice and regenerated through evacuation by a mechanic 
pump, it is likely that there is a minor oxygen contam- 
ination that may promote a chemical alteration to the 
Fe304 phase below 700 K. 

After the sample was cooled down from 780 K with a 
rate of 30 K/min, the magnetic transition around 700 K 
disappeared and the magnetization at room temperature 
was reduced to about 30% of the initial value, indicating 
that the Fe3S4 phase was completely decomposed into 
other phases (e^., Fe304 and nonmagnetic ferrosulfides 
such as FeS2 [8|) after the high-temperature magnetic 
measurement. 

In order to see the magnetic transition more clearly, 
we plot the derivative (dM/dT) of the warm-up magne- 
tization with respective to temperature. It is apparent 
that a minimum in dM/dT occurs at a temperature of 
about 677 K. The temperature corresponding to the min- 
imum in dM/dT or the inflection point of the M{T) curve 
should be the Curie temperature of greigite, that is, Tq — 
677 K. There is also a second local minimum at about 599 
K, which happens to the same as the Curie temperature 
of FcySs (Ref. nil). This suggests that a small fraction 
of the Fe3S4 phase was converted to the FeySs phase even 
at a temperature below 599 K. Since the oxygen contam- 
ination in our VSM system is so minor and the warming 
rate is high (30 K/min), the remaining Fe3S4 phase is still 
substantial when the Curie temperature of the phase is 
reached. This allows us to determine the Curie temper- 
ature of Fe3S4 quite reliably. The magnetic transition at 
677 K should be associated with the ferrimagnetic tran- 
sition of greigite unless Fe3S4 would have been converted 
to an unknown magnetic iron sulfide phase. 

Our data suggest that the intrinsic Curie temperature 
of greigite is about 677 K. The intrinsic net moment of 
greighte has been determined to be about 3.4 per 
formula unit from both magnetization [l^ and neutron 
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FIG. 3: a) Temperature dependence of the magnetization for 
sample A, which was measured in a field of 10 kOe. b) Tem- 
perature dependence of the derivative of the warm-up mag- 
netization. 



[14| data. Both the Curie temperature and the magnetic 
moment of greigite are significantly lower than those of 
its counterpart Fe304. 

In order to quantitatively understand the differences 
between the two spinels, we need to consistently extract 
the exchange energy Jab between the A and B sublat- 
tices for these two compounds. One way to determine 
Jab is from the temperature dependence of the satura- 
tion (or high-field) magnetization at low temperatures. 
Since exchange energies {J a A and Jbb) between ions in 
the same sublattice of the inverse spinel structure are 
negligibly small [IQi] , we can neglect these exchange ener- 
gies and only keep Jab- Within this approximation, the 
dispersion relation for the acoustic magnon to order fc^ 
is Q: 



Tiuj 



22 JabSaSb ,2 2 
16 2Sb - Sa ° ' 



(1) 



where w is the frequency and k is the wave-number of 
the spin wave. From this dispersion, one can readily ob- 
tain a relation for the saturation magnetization at low 



FIG. 4: a) Low-temperature saturation magnetization ver- 
sus T^/^ for a Fe304 crystal. The data are extracted from 
Ref. [20I ]. b) Low-temperature saturation magnetizations ver- 
sus T^/^ for two Fe3S4 samples (S706 and S504). The data 
are taken from Ref. 1121. 



temperatures: M,(T) = Af,(0)(l - CT^/^) with [gC 



C 



0.05864 1&{2Sb - SA)kB 



3/2 



4(25b — Sa) 22JabSaSb 



(2) 



where fcs is the Boltzman's constant. 

In Figure 4, we plot low-temperature saturation mag- 
netizations versus T'^/'^ for a Fe3 04 crystal (Fig. 4a) and 
for two Fe3S4 samples (Fig. 4b). The data for the Fe304 
crystal are extracted from Ref. I20II and the data for two 
Fe3S4 samples are from Ref. [1^. The best linear fits 
to the data yield C = (8.9±0.3) x lO'^ R-i '^ for Fe304; 
C = (1.15±0.06)xl0-5 for sample S706; C = 

(1.42±0.07)xl0"^ for sample S504. The average 

C value for Fe3S4 is calculated to be 1.29x10-^ K-^-^. 

For Fe304, 5"^ = 2.5 and Sb = 2.25 (Ref. [H). Sub- 
stituting C = (8.9±0.3)xl0-6 K-i •^ Sa = 2.5, and Sb 
= 2.25 into Eq. (2) yields JAB/kB = 22.7±0.6 K, which 
is slightly below the value (27.8 K) determined directly 
from neutron scattering [l9[. For Fe3S4, the two sub- 
lattice moments have also been determined by neutron 
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diffraction 14 1. From the measured values 



ruA 



3.08 fjLB and tob — 3.25 jiBi we find = 1.54 and Sb 
= 1.63. Substituting Sa = 1-54, S'b = 1.63, and the av- 
erage C value of 1.29x10"^ K~^-^ into Eq. (2), we find 
JAB/kB = 37.7 K. 

With neglect of exchange interactions between ions in 
the same sublattice of the inverse spinel structure, the 
Curie temperature within a mean-field approximation is 
given by [21| : 



T^^^ = 4V2{JAB/kB)VSASB{SA + 1){Sb + 1) (3) 

Since the mean-field theory tends to overestimate the 
Curie temperature. Corrections to the mean-field theory 
were given for different structures such as simple cubic 
(sc), body-centered cubic (bcc), and face-centered cubic 
(fee) j22|. Since the magnetic structure for the inverse 
spinel is fee [l§], the theoretical Curie temperature Tc 
after the correction to the mean- field theory is 22| : 



0.789[1 - 0.13/^SaSb{Sa + 1){Sb + 1)]tS'^. 



(4) 

For Fe304, Sa = 2.5, Sb = 2.25, and JAs/kB = 
22.7±0.6 K. Substituting these numbers into Eqs. (3) and 
(4), wc find that Tc = 808±22 K, which is slightly be- 
low the measured Tc = 851 K for Fe304 (Ref. 0). For 
Fe3S4, Sa = 1.54, Sb = 1.63, and JAs/ks = 37.7 K. 
Pluging these values into Eqs. (3) and (4) yields Tc — 
667 K, slightly below the measured value of 677 K. This 
quantitative agreement suggests that the sublattice mo- 
ments obtained from the neutron diffraction [l3| and the 
Curie temperature measured in this work represent the 
intrinsic magnetic properties of greigite. 



Now we compare the first principle calculations 15 1 
with the experimental results. Comparing the measured 
lattice constant a = 9.8538 A at 10 K (Ref. fl3|) with 
the calculated results [l^ yields U^/f = 1-16 eV, 
— 3.05 /is, and niB = 3.26 ^b- The predicted sub- 
lattice moments are in quantitative agreement with the 
measured values [3]: niA — 3.08±0.08 and niB — 
3.25±0.08 /iB. This quantitative agreement implies that 
the first principle calculations based on the GGA -I- U 
model is very reliable. This also implies that the inferred 
parameter Ueff = 1.16 eV is intrinsic to Fe3S4. This 
parameter Ueff = 1-16 eV is very important since it en- 
sures greigite to be a half metal [ll]. Compared with 
other half metals such as Cr02 (Ref. (2^) and doped 
manganites [24-26], greigite should be better for spin- 
tronic applications because it has a much higher Curie 
temperature than those of Cr02 and doped manganites. 
Compared with the half- metallic magnetite (Ref. [27|). 
Fe3S4 is also better for spintronic applications because it 
exhibits a highly conductive metallic behavior '^2§\ in con- 
trast to semiconductor-like electrical transport in Fe3 04 
[29| . Therefore, greigite is an excellent half- metallic ma- 
terial for spintronic applications. 



In summary, we have determined the intrinsic Curie 
temperature (677 K) and exchange energy (3.25 meV) 
between the two sublattices of greigite. The mesured 
Curie temperature is in quantitative agreement with that 
calculated using the exchange energy and the spin values 
of the two sublattices, which are independently deter- 
mined from the neutron and magnetization data of high- 
quality pure greigite samples 12, llJ] . We further show 
that, with an effective on-site Hubbard energy Ueff = 
1.16 eV, the lattice constant and two sublattice spins 
predicted from ah initio density-function theory are in 
nearly perfect agreement with the measured values. The 
parameter Ueff = 1-16 eV ensures Fe3S4 to be an ex- 
cellent half-metallic material for spintronic applications. 
The current work will benefit to scientists working in 
multiple disciplines including physics, materials science 
and technology, geophysics, geochemistry, and biology. 
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